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- . . In this paper we shall study whether dissipation in a X^p^ may be described, in the 

' long wavelength, low frequency limit, with a simple Ohmic term K(p, as it is usually 

, done, for example, in studies of defect formation in nonequilibrium phase transitions. 

O I We shall obtain an effective theory for the long wavelength modes through the coarse 
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graining of shorter wavelengths. We shall implement this coarse graining by iterating 
a Wilsonian renormalization group transformation, where infinitesimal momentum 



^ i shells are coarse-grained one at a time, on the influence action describing the dissi- 

' pative dynamics of the long wavelength modes. To the best of our knowledge, this is 



the first application of the nonequilibrium renormalization group to the calculation 
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, of a damping coefficient in quantum field theory. 
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Oh: 

^ I. INTRODUCTION 

;_i ■ In this paper we shall study whether dissipation in a X^p^ theory may be described, in the long 

^ wavelength, small frequency limit, with the addition of a simple Ohmic term KLp, as it is usually 
done, for example, in studies of defect formation in nonequilibrium phase transitions. We shall 
obtain an effective theory for the long wavelength modes through the coarse graining of shorter 
wavelengths. We shall implement this coarse graining by iterating a Wilsonian renormalization 
group (RG) transformation [H, 0, 0], where infinitesimal momentum shells are coarse-grained 
one at a time, of the influence action describing the dissipative dynamics of the long wavelength 
modes. To the best of our knowledge, this is the first application of the nonequilibrium RG to 
the calculation of a transport coefficient in quantum field theory. 

Understanding damping is one of the main goals of nonequilibrium quantum field theory, 
particularly with respect to (probably) its main applications, namely the generation of pri- 
mordial fluctuations during Inflation ^], the modelling of reheating after Inflation 5|, and the 
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formation of defects in nonequilibrium phase transitions [6[. There is a tension between the 
needs of model builders, who, understandably, seek simple and model-robust solutions, and first 
principles calculations, whose generic result is that damping in field theory is a complicated, 
non Markovian phenomenon 
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A clear depiction of this tension can be seen in the 
contrast between the early models of reheating, based on a linear damping of the infiaton field, 
and the now accepted picture of preheating, where the main decay mechanism for the infiaton 
is the parametric amplification of matter fields. This phenomenon is exponentially suppressed 
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in perturbative calculations 

The physical mechanism for damping in the long wavelength sector is the interaction with 
shorter wavelength modes. Damping is a feature of the effective theory where the shorter modes 



have been coarse-grained away 12j. Since this operation will leave the long wavelength modes 



in a mixed state, the natural description of the relevant sector is in terms of a density matrix 



131], and the natural action functional encoding the effective dynamics is the Feynman- Vernon 
Infiuence action iJ]. Now suppose we are given the infiuence action when all modes k > A 



have been coarse-grained away, and we wish to further coarse grain the modes in the range 
A > k > ko. A straightforward perturbative functional integration (for example, to two-loops 
accuracy) is likely to miss the most important effects. However, if we split the desired range 
into shells of infinitesimal thickness 6s, then each shell can be integrated out exactly within 
one-loop accuracy, because each loop integration brings a factor of 6s. Adding a change of units 
after each integration, we transform the shell coarse-graining into a RG fiow in the space of 
infiuence actions. Because we shall not assume equilibrium conditions, this may be called the 
nonequilibrium renormalization group. Our goal is to study the fiow of the damping parameter 
K in the infrared limit. 

The nonequilibrium renormalization group should not be confused with the so-called dynam- 
ical RG , which is concerned with the long term behavior of solutions to evolution equations 
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3,[3,[l8|, 



191]. The nonequilibrium RG, on the other hand, is closely related to the Martin, 



Siggia and Rose formalism 20|]. The equivalence between these two latter approaches has been 



discussed elsewhere 21, 22 



It is also important to stress the difference between our approach and the RG applied to 
dynamic phenomena, as in 23|, |2J, |25|, |26|, |27|, |28j . In these applications, a dissipative dynamics 
is assumed from scratch, as well as a phenomenological noise term to induce the fluctuations 
to be coarse grained away. In our application, instead, we shall assume that the initial point 
of the RG trajectory corresponds to a nondissipative theory. The source of noise will be just 
the quantum fluctuations in the short wavelength modes, as well as any statistical noise coming 
from the initial density matrix (see below). The issue is whether the renormalization group flow 
alone brings in dissipation and noise. 

It is important to stress two basic differences between the nonequilibrium and equilibrium 
renormalization groups 29|, ISOj. The influence action may be regarded as an action for a theory 
deflned on a "closed time path" (CTP) composed of two branches 3l(]. The flrst branch goes 
from the initial time t = to a later time t = T when the relevant observations will be performed; 
that is why we need the density matrix at T. The second branch returns from T to 0. Thus 
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each physical degree of freedom on the first branch acquires a twin on the second branch -we 
say the number of degrees of freedom is doubled. The influence action is not just a combination 
of the usual actions for each branch, but also admits direct couplings across the branches. The 
damping constant k is associated to one of these "mixed" terms. Therefore, the structure of the 
influence action (from now on, CTP action, to emphasize this feature) is much more complex 
than the usual Euclidean or "IN-OUT" action. 

The second fundamental difference is the presence of the parameter T itself. In nonequilibrium 
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evolution, it is important to specify the time scale over which we shall observe the system 
The CTP action contains this physical time scale T. From the point of view of the RG, this adds 
one more dimensional parameter to the theory, much as an external field in the Ising model. 
Physically, because time integrations are restricted to the interval [0,T], energy conservation 
does not hold at each vertex. This is of paramount importance regarding damping. 

The Wilsonian RG for equilibrium quantum fields in the imaginary time formalism was 
studied in Ref. 33], and in the real time formalism in Ref. j29|. The RG for the CTP effective 
action (obtained by taking the limit T oo) was studied by Dalvit and Mazzitelli 3J]; see also 



12] and 131] • Unlike those works, we focus on the dissipation and noise features of the effective 



dynamics, rather than in the running of the effective potential. For an application of similar 
ideas in a different field see 
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We shall now describe the basic assumptions underlying our work. 

In formulating a nonequilibrium RG, we must deal with the fact that the CTP action may 
have an arbitrary functional dependence on the fields and be nonlocal both in time and space. 
In principle, one can define an exact RG transformation [2I, Q, 36], where all three functional 



dependencies are left open. However, the resulting formalism is too complex to be of practical 
use 



3J, l37[. At the same time, one must beware of restricting the form of the action to the 
point of leaving out an important process. An example is, if we were to assume that the CTP 
action is the difference of an action functional for each branch of the CTP, we shall miss the 
damping and noise terms. This is allowed, of course, if one interest is, for example, the running 
of the effective potential, but it would be disastrous to our present concern. 

In our case, we will assume that the initial condition for the RG trajectory corresponds to a 
massless Ay?^ theory, with no noise or dissipation. This introduces a small parameter A in the 
model, and we will consider only the class of action functionals which may be generated from 
this initial condition to order A^. This means we shall include quadratic, quartic and six-field 
interactions. We shall allow the interaction terms to be nonlocal in time and in space, to the 
extend demanded by closure of the RG transformation. There are further constraints coming 
from the CTP boundary conditions; these will be described below. 

We shall work in three spatial dimensions. The structure of the nonequilibrium RG becomes 
much richer below this critical dimension, since there several fixed points (over and above the 
nontrivial fixed point of the equilibrium RG). These will be described in a separate publication. 
However, for the applications listed at the beginning, three spatial dimensions is the most 
relevant choice. We shall assume that the initial density matrix decomposes into an independent 
matrix for each mode. The formalism may be developed with great generality, but to obtain 
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definite results we will focus in the case were the initial density matrix corresponds to a free 
field at finite temperature. This is a nonequilibrium initial condition for the interacting field 
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We are mostly interested in the case of small T (that is, AT ^ 1), since for large T we expect 
the nonequilibrium RG to converge to the usual one (we shall return to this point below). This 
means that in computing Feynman graphs, it is allowed to use propagators pertaining to the 
initial density matrix and couplings, as any drift term will bring additional powers of A. To 
extend the nonequilibrium RG to a larger T range (for example, to study thermalization within 
this approach) a fully self-consistent approach is necessary 39|| in the spirit of the so-called 
"environmentally friendly" RG |40| . 

The paper is organized as follows. In the next two sections we provide an introduction to the 
CTP formalism, thus fixing our notation, and introduce the parametrization of the CTP action 
and of the initial density matrix. As we have remarked, the form of the CTP action is chosen 
to enforce closure of the RG transformation, within the constraints imposed by CTP boundary 
conditions. In Sec. [IV] we define the RG transformation. This is composed of two stages, first 
the elimination of one shell in momentum space, thus lowering the cutoff to (1 — 6s)A, and then 
the rescaling of momenta, times and fields to restore the value of the cutoff and the coefficients 
of the kinetic terms in the CTP action. We will show that, in spite of rescaling, the parameter T 
may be kept as a RG invariant. We then translate the change in the CTP action into a dynamic 
equation for the RG flow in parameter space. In Sec. |V]we write down these RG equations to 
order and discuss their structure. In Sec. IVTlwe discuss the flow of the damping parameter 
K and of the noise kernel z/ in the regime outlined above. 

We conclude the paper with some brief final remarks. We have concentrated most technical 
details into appendixes. 



II. OPEN SYSTEMS, CTP AND RG 

The core of the RG transformation is the elimination of selected degrees of freedom. When 
we are only concerned with the lower wave number sector of a field, we can carry out explicitly 
the integration over the higher wave number modes in the density matrix. As a result, these 
modes are eliminated from the description. This partial integration will return a description 
for the lower wave number modes only. The influence of the higher modes is incorporated in 
the parameters which define the effective action for the surviving modes. The procedure may 
be seen as a straightforward application of the Feynman- Vernon influence functional techniques 
for open systems. In the problem at hand, the long modes are regarded as the system and the 
short ones as the environment. 

We first review briefly the ideas behind the Feynman- Vernon formalism. We will deal with 
a system composed of two parts, which we call the relevant system, or simply the system, and 
the environment. 
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A. Open systems and Feynman- Vernon formalism 

Let if be the field variables of the whole system, and S its action. Then, its density matrix 
p admits the following representation 



p[<^+,<^-,T] = (<^+|f/(0,T)pf/(T,0)|<^-) 



V^+V^- e^Y>^{S[^^]-S[^-]) p[^+(0), ^-(0), 0]. (1) 

ip+{T)=v+,ip-{T)=ip- 

This expression says how to obtain the matrix elements of the density matrix operator in the 
Schrodinger representation at time T. The generating functional for the expectation values is 
obtained by adding current terms, and by closing the path of integration to take the trace 

Z[J+,J1 = jv^+V^~e^^ii^S[^+]-S[^-] + y"(JV+ + J-v^-)|p[v9+(0),^-(0),0]. (2) 

V+(T)=v-{T) 

Derivatives of Z with respect to the currents give expectation values of products of the fields 
(with a time ordering that depends on which derivatives have been taken). Equation ([2]) can be 
thought as an integral over single histories defined on a closed time path. This path has a first 
branch from to T, where the history takes the values ^^{t), and a second branch from T to 0, 
where the history takes the values f^{t). The CTP condition 

^HT)=^~{T) (3) 

implies that each history is a continuous function of the time along the path. The object 

ScTPy+,^-] = S[cp+]~S[^-] (4) 

is called the CTP action. In general it can be written 

Pb+,¥P-,T] = J Vip+Vif- exp{tScTp[^^,^-]) pb+(0), (/."(O), 0]. (5) 

Now, we make the division between relevant system and environment. Suppose that the 
whole system consists on two kind of variables and <^>, being the system variables, and 
(py the environment variables. Depending on the context, this could be a division between 
system and bath fields in the thermodynamical sense, or between different particle fields, or, as 
will be in our case, two different sectors of the same field: the lower and the higher wave number 
sectors. 

The whole system is described using the generating functional associated with the full density 
matrix, which depends on both system and environment variables. However, if we only wish 
to compute expectation values for the system observables, we can use the generating functional 
associated with the reduced density matrix that is obtained after taking the trace with respect to 
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the environment variables. To do this, suppose that the CTP action of the system + environment 
is ScTp[^p~^ ,^p~], then write it as a functional of the two pairs of fields, {ip^,{p~) and 
for the system and the environment, respectively, 

ScTp[^^, = 5'cTp[v5<, V5<] + ^S[(p%, V5<, V2>, (pZ.]. (6) 

For future convenience, the same functional S'ctp appears in both members; this is just the 
definition of AS*. The functional AS contains the information about the environment itself 
and about the interaction with the system. Let us further assume that at t = system and 
environment are uncorrelated, 

p(0) = p.(0)®pe(0), (7) 

where Ps(0) (pe(0)) refers to the initial density matrix of the proper system (environment). 
After, replacing ([6]) in the Eq. ([5]) and integrating out the environment variables ip^, we get an 
evolution law for the so called reduced density matrix 

Prb<,¥'<,T] = j V(p%Vip^ expi (^CTpbo + ^/Fb<,^<]) psi^ti*^), ^<{0),0], (8) 

¥J+(T)=<p+,</p-(T)=<p- 

where the influence action Sjp is given by 



^iSiF[<fit.''P<] = 



J Vip+Vip- exp^(A5[<^+,^-,<^+,yp;]) peb^(0),<^>(0),0]. (9) 



All the influence of the environment on the system is encoded into 5*/^. The elimination of the 
fields ipy takes the initial CTP action, S'ctp, and returns S'ctp<5 the coarse-grained CTP action 
for the proper system: 

5'cTP<b<,V5<] = S'cTpb<,V5<] + 5'/Fb<,V5<], (10) 

which is not necessarily the difference of some functional evaluated in each branch, S'[<y9^] — S'[y9^] 
for some S as in Eq. (jlj), but usually a more complex functional. It can entangle the two 
branches. CTP actions of this general type satisfy certain properties, namely 

S'cTpb,V5] = 0, (11) 

•S'ctp b",V5"^] = -^ctp^"^, V^"]*- (12) 
Thus, if we introduce new variables 



= - , (13) 

ip = p>+ + p>-, (14) 
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the expansion of 5'ctp in powers of 0, with functional coefficients (T„[(y9; ti, . . . , t„] depending on 
Lf and n time variables 

5'ctp [0,V?] = / rfti . . .(it„ a„[v9;ti, . . . ,t„] . . .0(t„), (15) 

n>l 

should start from = 1, and all the odd (even) terms should be real (imaginary). In what 
follows the fields and will be used instead of y?"*" and the CTP condition ([3]) is now 
expressed as 

0(T) = 0. (16) 

The application of this formalism to the case where the system and the environment are two 
sectors of the same scalar field is presented in Appendix [XI 



III. PARAMETRIZATION OF THE CTP ACTION AND OF THE INITIAL 

DENSITY MATRIX. 



In this section we introduce the parametrization of the action and of the initial density 
matrix. The parametrization of the action will be quite general. In Sec. IIVI we will define the 
RG transformation for this general class of actions, without any further assumption regarding the 
initial condition at cutoff A. In Sec. |V]the initial condition will be fixed and an approximation 
scheme introduced. This will allow us to reduce the parameters to be considered to a finite 
number of couplings with simpler dependencies on the momenta and time variables. 

We will use natural units in which A = 1. To compute the RG transformation we will assume 
that the action is 



S'[0, if] = So[(j), <f] + S'int[0, (/?], 



(17) 



where the free action is 



dt / d\ 





i0(k, t) vi(-k, t) - ^0(k, t) {e + m^) <^(-k, t) 



and the interaction part 



For n = 2 



-K 0(k,t)0(-k,t) + -V 0(k,t)0(-k,t) 



even n >2 



[ dti I dt2 I d% Li(fc;ti,t2) 0(k,ti) vp(-k,t2) 

JQ JO J A '- 



+ i V22{k; ti, t2) 0(k, ti) 0(-k, t2) 



(19) 
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and for n > 2 



1—- 




X 



-(-l)-]/2 



({k};{t}) 



(20) 



m=l 



The factor fi^ where is the area of the unit sphere in d dimensions, is introduced for 
further convenience. The subscript A in the integral symbols means that the integrations are 
restricted to IkA < A = 1. Moreover, we have defined 



({k};W) '^nm(kl, • • • , k^, . . . , k^, 

ti, . . . ,tm, ■ ■ ■ ,tn) (with l<m<n), 

^n-m ^ ^(]^^^ ^^y _ _ ^(j^^^ (^(k„+i, Wl) . . . ip{K, Q . (21) 



With these definitions, both the fields and the couplings are dimensionless. Note that the CTP 
action associated with the usual classical action for a massless X(f)'^ is 



/ dt 

48 70 



dt / d'^k - 



0(k,t) 0(-k,t) - A;V(k,t)</?(-k,t) 



d% . . . d% 



(ki, t)v3(k2, t)v5(k3, t)v5(k4, t) 



. 1=1 



+0(ki,t)0(k2,t)0(k3,t)(/^(k4,t) 

Hence, in the notation just introduced, only v^i and ^43 are different from zero, 



V41 = V43 



48(27r)' 



A5(ti-t2)5(ti-t3)5(ii-t4). 



(22) 



(23) 



(The factor A'^"^, actually equal to 1, has been included for completeness.) 

According to the definitions (fT9|) and (!20|) . Vnm couples n fields altogether, m of type and 
n — m of type if, always with m > 1. In addition, when in a given term the number m of 
fields of the type (p is even, the coupling function appears multiplied by i. Therefore all the 
couplings functions Vnm will be real [see Eq. f|T2|) ]. We will take the couplings symmetrical with 
respect to the permutations of the variables of each type of field; e.g., U42 is unchanged when 
(ki,i:i) (k2,i:2) or (k3,t3) (k4,t4). We will also assume that the Vnm couplings do not 
contain time derivatives of Dirac deltas. 

To make sure that the RG transformation relates actions on the same class and preserves the 
CTP properties ( fTTI) and ( fT2l) . the couplings Vnm must be such that at least one of the fields 
in each interaction term is evaluated at a time t equal or later than the fields (f. This condition 
is directly related to the fact that the expectation value {ip(t)(j)(t')) is causal. 
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Finally, we give a prescription to separate out the free action from the quadratic interaction 
terms. It must be 



[ dt [ dt' V2ii0; t,t') = 0, 
Jo Jo 



Q2 / rT .T 

— (I dt I dt' V2iik;t,t') 



dk 



JO 



k=0 



dt j dt' {t-t') V2i{k]t,t') = 0, 
Jo 

[ dt [ dt' V22{k;t,t') = 0. 
Jo Jo 



(24) 



In this way, the terms m'^(j)ip, k'^cpip, Kcfxp and vcjP' are isolated within Sq. We will return to this 
point in Appendix O 

We have to chose the initial conditions p[0(k, 0), V5(k, 0)]. We will use 



p[0(k,O),<^(k,O)]=exp<^ - / d\ 



tanhl^^^^M^ ^fk,0)(/.(-k,0) 



+ coth(«) 



4 \ 2 ) 

^(k,0)<^(-k,0) 



(25) 



This corresponds to uncoupled, free fields with frequency a{k) and with a temperature [3{k)~^ 
which depends on the wave number. Observe, however, that this is a nonequilibrium initial 
condition for the interacting fields. 

The propagators corresponding to the free action fllSp and the density matrix (123]) are given 
in Appendix [Bl 



IV. THE RG TRANSFORMATION 

In this section, we define the RG transformation for the general class of actions f|T7|) . The 
infinitesimal RG transformation, from which the finite transformation is obtained, is composed 
by two operations: i) elimination of the modes with wave numbers in an infinitesimal momentum 
shell, and ii) rescaling. We analyze each of these operations in the following two subsections. 
The joint effect of i) and ii) is written as a set of differential equations in Sec. IIV CI 



A. Mode elimination 



The first step to obtain the RG equations is discussed in Appendix \^ Here we will apply 
the results of that section to the action f[T7|) . Modes with momenta in the shell between h and 
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1 are eliminated {b = 1 — Ss, < Ss 1). After eliminating the modes in the shell, one obtains 
an effective action for the modes with momenta below b. The effective action is given by Eqs. 
[9| and flAlOl) . It can be written as 



^< = ^0 + SL (26) 

where 



T 

JbA 



-{k + 5k) (j)ip+'^-{v + 5v) (t)(t)), (27) 
and where 5*;^^ has the same form that 5'int in Eq. (fT8|) . but with perturbed couplings 

'"'nm = ^nm + 5Vnm- (28) 

Note that the momenta in the integrals now go up to bA. We discuss in the Appendix O the 
actual calculations involved. 

At this stage, the initial density matrix remains unchanged. 

B. Rescaling 

Now we introduce the second step in the RG transformation. We redefine the fields and 
change integration variables in the action, writing 

b"^(j){b-%b'^H) (29) 

and 

(&"^k, (30) 

instead of 0(k, t) and ip(k,t), and 

6k (31) 

and 

b-^H (32) 

instead of k and t. The rescaling redefines fields, times and momenta and restores the cutoff to 
its original value A = 1. The two exponents, 

= ^{r] -d-1), (33) 
at = l-v, (34) 
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are chosen to reduce to 1 the coefficients of the terms (fxp and fc^c/x/? in the free action. Here 77 
is the quantity introduced in fl27j) . Rescahng affects both the couphngs in the action and the 
parameters which define the density matrix. 

It is important to notice that, in spite of rescahng, the time T may be kept as a RG invariant. 
This is a consequence of the CTP condition ( fT6ll and causahty. This fact and other details 
concerning rescahng are shown in Appendix [Dl 



C. RG equations and formal solutions 



The effect of mode ehmination and rescahng over the parameters which define the action and 
the density matrix can be summarized in the following differential equations 



— - 2 + 2r/ = -—, 

OS J OS 



(35) 



1 k ^ ^ - 

ds ^ dk ) ^ 6s' 



(36) 



d 2 + 2+A;^^ 

ds dk J 6s' 



(37) 



d /n \ , „ 3n 



El , d d 
1 ^'dk, - 



i=l 



6Vr, 



6s 



{31 



d , , d \ ^ 



0. 



(39) 
(40) 



[See Eqs. flD4l) - flD10p . Equations. fl33l) and fl34l) have been used for the exponents.] 

Equations ( !35|) -(l40l). together with the equation for r], are the RG equations. Integrating 
these differential equations is equivalent to iterate the infinitesimal transformation. They have 
the general form 



d_ 

ds 



+ ap + ^ ( ki 

i=i ^ 



— - t — 

dki * dti 



(41) 



where ap and at are functions of s alone. The formal solution F({k}; {t},s), with the initial 
condition 



F({k};{t};0)=Fo({k};{t}), 



(42) 
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IS 



F({k}; {t}, s) = J' ds' e^^(^')-^^(^) g (^e''"{k}; e^'^^)-^*^^'^*}, s') + 
e-/3FWFo (e-^{k};e^*(^){t}) , 



(43) 



where Pf{s) = / ds' ap{s'), and = / (is' Q;t(s'). 
Jo Jo 



V. THE RG EQUATIONS TO ORDER 

A. Reduced set of parameters 

So far, the RG transformation defined for the class of actions ( fTTIl is rather general. The 
transformation is closed with respect to this class, but infinitely many couplings have to be taken 
into account. To reduce the number of couplings that have to be considered, we will choose as 
the initial condition at the cutoff A the action ( |22l) . and compute the RG equations to order 
A^. The central question is to find the minimum set of parameters such that the transformation 
is closed to order A^ for this initial condition. This set, which has 18 elements, is found in 
Appendix [El together with a number of constraints upon the dependencies on the momenta and 
time variables. 

These results can be summarized by saying that, starting at s = from the action fl22|) . to 
order A^, the action for s > can be written as S* = 5*0 + ^int, where 5*0 is the free action of Eq. 
(fT8l) . and where Si^t is given by 



^int= [ dt [ d%V2l{t,s) (f){k,t)ip{-k,t) 
Jo J A 

dt dt' j d% ^W2i{k;t,t',s) (p{k,t)ip{~k,t') +iW22{k;t,t',s) 0(k, t)0(-k, t') 

I dt I Wd%5'^{ki + --- + ki) \Vii{s) 01<^2<^3<^4 + V^43(s) 010203<^4 



+ 



i=l 



/ dt / dt' / \{d%5''{ki + --- + ki) Wii{Qi2\t,t',s) 01(^2 vWi 

Jo Jo J A '- 

+iW42iQ 12; 't,t',s) ^1^2 (p'3'f'4 + W^3iQl2;t,t',s) 0iV?2 (^'^(p'i 

dt' / J]^ (ilfci (ki H hke) V6i{Qi23;t,t' , s) 0iV22V53 ^Wb^'e 



f / dt 
Jo JO 



i=l 



+'i'V62iQ 123; 't,t',s) 01V32V53 (l)'4 f 5^6 + '^S (Q 123; t^t'^s) 010203 ^Wb'^G 
(Q 123] t,t',s) (f)i^2V3 (PWb^G + ^^Q^iQ 123] t,t' , s) (f)i^2^3 0405^6 
+V65{Ql23; 't,t',s) 010203 0405V56 + *^66 (Ql23 ; i^, t', s) 010203 040506 



(44) 
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Here 0, = 0(kj,t), = 0(kj,t') (same for ip), and 

Qi2... = |ki + k2 + ...|. (45) 

Moreover, m^, V21, V41 and V43 are of order A, while k, z/, the rest of the couphngs in Sint, and 
rj are of order A^. 

It is straightforward to recast the general analysis of Sec. IIVI for this particular action. The 
RG group equations are given in the next section. 



B. The RG equations 

In the previous subsection we gave the reduced set of parameters that has to be considered 
in working out the equations to order A^. Now the RG equations will be grouped in the same 
order that will be solved later on. The groups form a hierarchy. The equations at the top are 
themselves a closed system. Once their solutions are known, the equations at the following level 
can be considered as effectively closed, and so on. 

At the top of the hierarchy are the equations that determine the propagators Q and G. It is 
important to notice that in writing the RG equations to order A^, for most cases it is enough to 
know the propagators Q and G to order zero in A. At that order they are independent of the 
couplings and are determined by the functions a and (3 in Eq. ( l25l) alone. The only exceptions 
are the equations for V21 and m^, where Q must be known to order A, in which case it will 
depend on m^. But these equations are at the end of the hierarchy, so they do not affect the 
equations that precede them and still will form a closed system. 

Thus, for the moment, we only need the propagators to order zero in A evaluated at A; = A = 1. 
From their general expressions ( IBll) and ( lB2p we find 

G{t, t') = 2 sin(t - t') e{t - t'), (46) 

and 

g{t, t', s) = — ^ [1 + 2/ (a (3{1, s)]] {a(l, sf cos(t - t') + \l - a(l, s)^] cos(t + t')}. (47) 
a(l, s) L V / J ^ 

The functions a(l, s) and /5(1, s) to order zero in A are given by Eqs. (139|) and (HO!) setting 77 = 0. 
At the second level are the equations for V41 and V43 

Ys+d-^^V,,is) = 18V,,isf{Gg}is,T), (48) 
^ + d-3]v43{s) = 18V^iV43{s){Gg}{s,T), (49) 



ds 



where 



{Gg} = viGg). (50) 
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It has been used that VQi{l;t,t' , s) = for the chosen initial conditions, and that t] can be set 
equal to zero in the left-hand side member of the RG equations to order [see Eqs. (|35ll-(|38l)]. 
Note that if V4i(s) = V43(s) for s = 0, then this relation holds for s > as well. The solutions 
V41 and V43 will depend on T through the function {G^}. 

At the next level we find the equations for the Wgi couplings, Appendix [El They have the 
same general form; we write only the first three equations, that are all that is needed for present 
discussion 



{'D + 2d-5)v6iik;t,t',s) = 3Vii{sf G{t,t') 6{k - V 



(51) 



(D + 2rf - 5) ve2{k; t, t', s) = ^V^iisf g{t, t\ s) 6{k - 1+), 



(52) 



(D + 2c/ - 5) v^^^{k] t, t', s) = 9V4i{s)V43{s) G(t, t') 6{k - 1+). 
[See Eq. (]D14p for a definition of 5 (A; — 1+).] Here 



^ d , d f d , d , 



(53) 



(54) 



According to the results of Appendix [Ej the factor at, which would have to appear in front of 
the time derivatives, has been set equal to one. 

Next come the equations for W41, W42 and W43, which are 



(D + ci - 4) Wiiik; t, t\ s) = 184;o V^i^sf GQit, t', s) - 18^4^ {G^;} 5{t - t') 
dfld 



+ 



fin 



[Gv(ii{Qkn; t, t', s) Q{t, t', s) + 2vti2{Qkn; t, t', s) G(t, t')] 



(55) 



D + d - 4) 1^42 (A:; t, t', s) = 95,,o V^4i(s) V^4i(s) GGit, t', s) - 2^/43(5) [GG(t, t')]^ 



a7 



v^KQkn] t, t', s) G{t, f) _ + AvMkn; t, t', s) g{t, t', s) 



(56) 



[B + d-A) W^sik; t, t', s) = l8Sk;o VM^is) GQit, t', s) - ISV^^V^sis) {G6?}(s) Sit - t') 



+ 



dVta 



2vQ^{QkQ; t, t', s) Q{t, t', s) + QvQiiQkQ] t, t', s) G{t, t') 



We defined 



Q 



\k + Vt\. 



(57) 



(5^ 



The subscript W means symmetrization respect to t and t' (we used that Vq2 and Q are already 
symmetric in these variables). The terms in Eqs. (I55p and (IFT]) proportional to 5{t — t') are the 
necessary subtractions after isolating the contributions to dVn and (5V43. 
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Finally, we write the equations for the the 2 field couplings and t], which can be solved in 
closed form if the solutions of the previous equations are known. We define three auxiliary 
functions 

fi{t,s)= [ dt'\w4i{0-t,t',s)g{t',t',s)+V4i{s)V2i{t',s)Gg{t,t',s)] +3V4i{s)g{t,t,s),{59) 
Jo ^ ^ 



[w^iiQkn; t, t', s) Git, t\ s) + W^42(QfcQ; t, t\ s) G(t, t') 



(60) 



[WMkn; t, t', s) g{t, t', s) + 2WMkn] t, t', s) G(t, t')] . (61) 



Each function is the sum of several diagrams appearing in Appendix [El 
In terms of these functions, it results 



{B-2)V2i{t,s) = h{t,s)-Vh{s) 



(62) 



(D - 3) 1^21 {k; t, t', s) = f2{k- 1, t', s) - 



5{t - 1') 



■Qf2ik,s) 2 



d 



^ + m - m 



Sit - 1') 



(63) 



(D - 3) W22ik; t, s) = h{k- 1, t\ s) - Vfsik, s) 6(1 - t'), 



(D - 2) m\s) = -2 [Phis) + Vf2{0, s)] 



(64) 
(65) 



(D-l)K{k,s) = -Qf2{k,s), (66) 
{B~2)u{k,s) = 2VMk,s), (67) 

.= -2^(0..). (68) 

The operators V and Q are defined in Appendix [Ul Eqs. (1C9P and (IClOp . Note also that, as 
was mentioned at the beginning of this subsection, since V41 is of order A, in the last term of 
Eq. ( l59i) . Q must be written to order A. In this way, the equations for V21 and will include 
all the terms to order A^. As we will not attempt to find the solutions for V21 and to order 
A^ we will not need to write Q to order A explicitly. 

The way to integrate the RG equations starts by writing down the propagators G and Q as 
functions of s. Next the equations for V41 and V43 can be integrated immediately. Once these 
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functions are known, the equations for the fgi couphngs are solved using the formal solution 
(H3l) . With these solutions, the equations for the can be integrated, and then the equations 
for V21, W2i, m? , n, V and r]. This is immediate, except for V21 and m^, because they are given 
by integro-differential equations. (They become ordinary differential equations to order A.) 

We will work in c? = 3. First we find the propagators to order zero in A. We need initial 
conditions for a and /3 in the density matrix. We choose 



a(fc, 0) = uj^ik, 0) = /c. 



and 



/3(A;,0) = /?o(A^). 

For the moment /5o can remain unspecified. From Eq. fl39l) . with 77 = 

a(A;, s) = k. 



(69) 



(70) 



(71) 



Note that the evolution of the product a(3, that appears in the argument of the function / in 
Eq. (]A17|) for Q, is given by 



which has the immediate solution 



ap{k,s) = ap{ke-',0). 



(72) 



(73) 



It is convenient to use the variable z = e'^ instead of s in the propagators. Using (|69|) and (|70 
in fl75]) . and then replacing in fll7|) . setting k = 1 



g{t,t',z)=2 l + 2flz'^po{z-^)) cos(t-t'), 



from which Eq. (|50|) gives 

{Gg}{z,T) 



1 sin(2T) 

2T ^ ' 



(74) 



(75) 



The dependence of {G^} on T comes through the first term in Eq. ( |75l) . In Fig. 1 we plot 
1 — (2T)~^ sin(2T) as a function of x = 1/T. Observe that for T 00 the function tends to 1. 
Then, from Eqs. fj^Hj) and (jUj), to order A^, we get 



V,,{s) = VUs) = Vo + 18^0^ ['ds' {Gg}{e'',T), 

Jo 



where 



{27rY 48 



d=3 



1 A 

'2^ 48' 



(76) 



(77) 
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FIG. 1: The graph of 1 - {2T)-^ sm{2T) as a function of x = T-^. The term {G^} depends on T through this 
function and has a definite hmit when T ^ oo. 



To solve the equations for the fgi couphngs to order A^, in the second members of Eqs. 
( l5T|) - fl53|l we can use the initial values of V41 and V43. Using fH3l) 



(78) 

(79) 
(80) 



vei{k;t,t',s) = -Vq^ G{k{t,t')) e(lV - fc)e(/t - 1+), 
ve2ik; t, t\ s) = ^ Q (k{t, t'), j) e(lV - k)Qik - 1+), 



,(1) 

^63 



9 



ik;t,t',s) = -rVo G{k{t,t')) Q{l+e' -k)Q{k-l+). 

rC 



[See Eq. f lESP for a definition of Q{k — 1"^).] The notation ( k(t, t') ) means (H, kt'), and will be 
used frequently. 

Next the equations for the are integrated: 



Wii{k-t,t',s) = 18 4;oK)'2:i[Ge?](t,t',e'^ 



+9^2 



2X2[G, G] {k- 1, t', e') + J2[G, g] (k, t, t', 



V,i{s)-V, 



^1) 



W^2{k;t,t\s) = 9 5k;oVo^ 



+ 18 



MG, Q\ {k- 1, t', e') - 2 J2[G, G] (fc; t, t', e^)^ 



(82) 
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+18 



J2 [e?, G] (A;; t, t', e^) + X2 [G, e?] (fc; t, t', e" 



where 



S{t-t'), (83) 



and 



c/fig ^2 (\z~^e''k + Q\ze-''{t,t'),\z~^e''k + Q\-^e'' 
Here, the angular integration is over the region S defined by 



1 < 



z-'^e''k + n 



^5) 



^6) 



so X2 = for /c = (though its hmit when A; — > 0"*" is not necessarily 0). As before, z is identified 
with e^. 

The function I2 can be expanded as 



l2[Gi,G2] = Q{l + k-z) X2<[Gi,G2] + e(^-l-A;) J2>[Gi,G2] 



^7) 



where 



12< [Gi, G2] = ^[dy f dx Gi (^y (t, t'), '-^ G2 (^x (t, t'), ^) 



and 



z—k ry+k 



X2>[Gi,G2] = ^|^y^~ dy j""^dx + j"^ ^dy j" dx^ Gi (^y (t, t'), ^ j G2 (^x (t, t'), • (89) 

X2 is continuous and have continuous first partial derivatives along the line z = k + 1. 

We end this section taking the limit of T — ^> 00, and zero initial temperature, i.e. Pq'^ = 0. 
We focus in the equations for A and m^. In this limit, Eq. (!75l) gives 



{G6;} = 1 (T^oo, f3^' = 0). 



(90) 



With initial conditions 



K4l(0) = 1^43(0) 



A 



(91) 
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the equations for V41 and V43 are 



d \ 
— + d-3] Viiis: 



(92) 



For s > we can define 



A(.) 



48(27r)'^ 

n7~ 



wfiich will satisfies the following differential equation 



d 



+ d-3]X 



ds J 8{2ttY 

On the other hand, to order A, the equation ( l65l) for is 



(1-^1 



A. 



In computing Vf\ to order A in Eq. fl65|) . it has been used that 

Pg = 2 (T^oo, /3o"' = 0). 



(93) 



(94) 



(95) 



(96) 



(97) 



Equations ( l95l) and ( !97l) are the CTP equivalents of the textbook equations derived using the 
RG group defined for an Euclidean action 4l|] . Observe that we are integrating out momentum 
shells of infinitesimal width in the spatial directions but leaving the time direction unrestricted. 
Thus the numerical coefficients in Eqs. (1951) and (1971) may be different from those obtained when 
a spherical shell in Euclidean momentum space is integrated out [33 . 



VI. THE DAMPING CONSTANT n AND THE NOISE KERNEL f 



A. The damping constant 

To compute k we need the function /2, Eq. (IBO]) . and then Q/2. Not all the terms in 14^41 
and W42, given in Eqs. flHTl) and fl82|) . contribute to k. The terms proportional to 6{t — t'), for 
which the application of Q gives zero, can be discarded, as well as the terms proportional to 
the Kronecker deltas, because they do not contribute to the angular integrals. If Eq. (l66l) is 
rewritten retaining only the relevant terms it yields 



(D - 1)k = -I8V0' Q/2, 



(98) 
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where 



f;ik;t,t',z) 



2l2[g, G] {Qkn; t, t',z)+ 12 [G, Q] {Qkn, t, t', z) 



g{t,t',z) 



I2 [Q: G] {Qkn; t, t', z) - J2 [G, G] (Q^o; t, t') 



G{t,t'). 



(99) 



In the last hne, it has been used that for the function G given in the Eq. (1461) is 
J2[G, G](fc; t, t', e^)H7 G(t, t') = ^I^IG, G](fc; t, t', e^) G(t, t'). 



(100) 



(101) 



Using the formal solution in Eq. (I43l) 

/€(fc, s) = -18V^ [ ds' e'-'' Qf;{e''-' k, e''). 
Jo 

Due to the general structure of the functions X2, Eq. (187|) . /t will be given by an expression of 
the form 



with 



+e (^e'' - I k e^'"^ + fi| - 1) A>(| k e''-' + e^') 



2J2< [g, G] + J2< [G, ^]) ^ + 2 (J2< - J2< [G, G]) G 



(102) 



(103) 



and with a similar expression for A>. Observe that both and A> depend on T through the 
application of Q, and that this dependence will be inherited by k. 
Changing variables, Eq. fll02p can be rewritten as 



2 r^+^ 



K(k, z) = 9Vn -r dy du ^ 



Q{u + l-y)A^ {u,y) 
+eiy-u-l)A> {u,y) 



(104) 



In doing this integral, there are three different regions of the plane kz, restricted to < A; < 1 
and 1 < z, that must be analyzed separately: i) 2; < 2 — fc, ii) 2 — /c < 2; < 2 + A;, and iii) 
2 + k < z. In each region k is given by the following expressions: 
Region i: 



K{k, z) = 9Vq — I dy I ^ du — A< {u, y) . 



1+^ 



u 



(105) 
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Region ii: 



2z 
z+fe 



dy 



1+ 



yk 



1- 



yk 



k 



du + 



2z 
z+k 



2z 
z + k 



u 



dy 



yk 



U 



rfw — v4< (u, y) 



+9K)V / / du^A>{u,y). 



^_yk y- 



(106) 



Region iii: 



K{k, z) = 9V^ 



2z 
z + k 



dy 



yk 



du + 



2z 
z — k 



dy f' 



z + k ^ 



yk 



du 



u 



2z 

z — k 



2z 
z + k 



dy 



du + 



yk 



2z 
z-k 



dy 



yk 



U 



yk 



du] ^ Ay {u,y) . 



(107) 



Note that for k = the angular integral in Eq. (11021) is trivial 

niO,z) = 18V,'z r^A<(l,y) 9(2 

Ji y 

r /-s 

+18V^ z 



1 



2 y 



%Ay{l,y) 



0(^-2), 



(108) 



which is also the limit when A; ^ of the expressions (11051) and (11071) for 2; < 2 and 2; > 2, 
respectively. 

For z^ 1 n becomes independent of k. If A^il^y) grows slower than y, for z 3> 1 it results 



k{z) ~ 1%V^ z 



' f % A.(l,y) + r % A,(l,y) 

Ji y^ J2 y 



Ko Z. 



(109) 



This expression corresponds to what can be expected from dimensional arguments. The constant 
kq can be thought as the k accumulated during the first stages of evolution. The factor z comes 
solely from the rescaling. If A>(1, y) grows as y or faster for y 00, to leading order, k will be 
given by the last term in Eq. (11081) 



k{z) ~ 181^0^ z 



2 y 



^^>(i,y). 



(110) 



Thus, in this case k, will acquire an anomalous dimension. 



B. The noise kernel v 

The noise kernel u is obtained following the same steps used in the calculation of k. First, 
Eq. (1671) is rewritten keeping in only the relevant terms. Term proportional to 6{t — t') in 
the solution (l83l) for W43 drops out because W43 appears multiplied by G{t,t'), which is zero 
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in t = t'). The terms proportional to 5k-fl also drop out because they do not contribute to the 
angular average. The equation for z/ reduces to 



where 



f*{k-t,t\z) = 2 



{■D-2)iy = 3QV,'Pf;, 



UG^ G]iQkn; t, t\ z) + J2[G, g]{Qkn, t, t\ z) 



(111) 



G(t,t',z) 



+ 



12[Q, Q]{Qm\ t, t', z) - J2[G, G]{Qkn\ t, t') 



(112) 



Defining 



B^=V 



2 (x2< [g, G] + X2< [G, e;]) G + (j2< [g,g]- 12< [g, G])g 



(113) 



and analogously i?>, it results 



u{k,z) = 18Vij 



dy 



du 



u 



y 

+e(y-n-l)5> {u,y) 



e{u + l-y) 5< (m,i/) 



(114) 



Again, note that the functions B will depend on T through the application of V, an so will v. 
Expressions analogous to Eqs. (llU5l) - (llU7p hold for v. For /c = 



v{Q,z) = ?,QV^z' r^B^{l,y) 6(2 

Ji r 



+36^2 z^ 



^dy ^ . f dy 



Ui y 



^B4l,y)+ / ^i?>(l,y) 



2 y 



Q{z-2). 



(115) 



For z ^ 1, V becomes independent of k. If B^{l,y) grows slower than y"^, the behavior of u is 
controlled by the rescaling 



u{z) ~ 36V^ z^ 



Ji y^ J2 r 

Otherwise i' acquires an anomalous dimension 

uiz) 36V,' z' r%B^{l,y). 

j2 y 



y) 



(116) 



;ii7) 
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C. K and V for a specific choice of (3q 



We compute k and u choosing /3(^^(A;) = 0, namely the case of zero initial temperature. The 
function 1 + 2/ in Eq. (17i|) . with / defined in Eq. (IBSp . is replaced by 1. We compute k and v 
only for = 0, since the dependence on /c for 2; 3> 1 is weak. The detailed expressions are given 
in Appendix [Fl 

The actual calculation shows that Ay(l,y) and B^{l,y) both grow as y when y ^ 00. This 
implies that, for 2; 3> 1, k will develop an anomalous dimension, while u will go simply as z"^. 
Making explicit the dependence on T, for z ^ 1 



k{z,T) 



9Vl 
AT 



zlogz [7 - 2T^ - 8 cosT + cos{2T)] 



and 



z/(fc,T) 
+ 4 InT- 16 



9V' 



135 + 4 



8T 

Stt + 3T - 16 Si(2T) 



7i5 - 34 In 2 - 7 In 3 + 7 Ci(3T) - 3 cos(4T) 



77r + 6T - 14 Si(4T) sin(2T) 



Ci(4T) - 8 cos(3T) 



sinT + 4 

- 12 19 - AT^ Ci(T) +8 7 - 6T^ Ci(2T) +4 1 - 4T^ 

- 24[7-8Ci(2T)] cosT + 4 [l 1 + 2 Ci(T) - 2 Ci(2T) - 14Ci(4T) cos(2T) 

+ AT [iSrr - 54 Si(T) + 28 Si(2T) - 6 Si(3T) - 2 Si(4T) + sin(4T)] |, (118) 

where Si and Ci are the sin and the cos integral functions. Observe that the damping constant k 
is not definite positive. This suggests that the underlying mechanism could be similar to Landau 
damping 421 . 



VII. FINAL REMARKS 

In this paper we have investigated the nonequilibrium dynamics of the low frequency, long 
wavelength modes of a self-interacting real scalar field. We have computed the influence func- 
tional encoding the back-reaction on these modes of the higher frequency, shorter wavelength 
sector. We have obtained the coefficients of the influence functional by solving the RG equations 
which describe the change in the influence functional induced by the progressive averaging over 
momentum shells. 

The main finding of this paper is that the infiuence functional for the low frequency modes 
contains terms associated with damping and noise. It is most important to stress that we 
have not put these terms by hand; contrariwise, we have assumed that they are absent at high 
energies. Damping and noise are forced upon us by the RG fiow itself. 

We have been able to retrieve these terms because we have gone beyond the adiabatic ap- 
proximation for the environment (high frequency) modes. A crucial step in this direction is the 
recognition of the role of the parameter T in nonequilibrium evolution. Because time-integration 
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is restricted to the lapse from to T, energy fluctuations, and thus particle creation from the 
vacuum in the environment sector, are allowed even in simple approximations to the dynamics. 
Although the flnite temperature RG in the real time formulation is well known, and also the 
running of the effective potential in a OTP framework, to the best of our knowledge this is 
the flrst instance were the RG is used to compute an intrinsically nonequilibrium feature in 
relativistic fleld theory. 

The most limiting approximation we have made is the use of flxed propagators in internal 
lines, disregarding the changes in the propagators caused (mainly, but not only) by the generation 
of mass, dissipation and noise terms through the RG flow. This approximation forces us to 
restrict our analysis to very short times (a few inverse cutoffs). A fully self-consistent RG should 
overcome this shortcoming. Also we have only considered a very simple -though nonequilibrium- 
initial condition for the fleld. A more flexible approach regarding initial conditions is required 
for most interesting applications. 

The methods advanced in this paper are relevant to essentially all the applications of nonequi- 
librium quantum fleld theory, from high energy ones like thermalization in relativistic heavy ion 
collisions and reheating after inflation to low energy applications such as the dynamics of glasses 
and Bose-Einstein condensates. It is clear that we have the bare outlines of a framework, but 
we look forward to see this framework develop and fructify in manifold ways. 
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APPENDIX A: INTEGRATING OUT MODES IN A SHELL 

In the problem at hand, the relevant system and the environment are sectors of the same 
scalar fleld, being 0> and the variables associated with the short scales. The elimination of 
these modes is the fundamental step in the RG transformation, and follows the lines sketched 
in Sec. UTAl 

The CTP generating functional is 

Z[3,j]=jv^V(P exp^ l^cTP [^,^]+j^dt jd'^k [J(k, t)(^(k, t) + j(k, t)0(k, t)]| p[0o, <^o]- (Al) 

Here the integral is over histories with 0(k, T) = 0, and (po and ipQ stand for the flelds at 
t = 0. We retain explicitly the dimensional parameter T at which both histories on the CTP 
integral coincide. The fleld conflgurations contain modes up to momentum A: 0(k, t) and fCk, t) 
are deflned in the region |k| = < A, and so are the currents J and j. If necessary, flelds and 
currents can be deflned to be zero outside this domain. The subscript in the momentum integral 
indicates that the integration domain is bounded by A; = A. 
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The density matrix is known at t = 0. For simplicity we assume that 

with the same functional p everywhere [cf. Eq. fl25] . 

Only modes with wave numbers within an infinitesimal shell will be integrated out. Let 
< 5s ^ 1 the infinitesimal parameter of the transformation, and define 

h=l-5s. (A3) 

The momentum domain /c < A is divided in two regions: the shell 

6A < A; < A, (A4) 

and its interior: 

k<hK. (A5) 

Accordingly, the fields are split in two parts, one which contains the modes within the shell, 
indicated by (/)>, and another part containing the modes with k < bA, indicated by 0<, i.e. 

= 0< + 0>, 

where 

0<(k,t) = 0(k,t)0(6A- A;), (A6) 
0>(k, t) = 0(k, t) [e{k - bA) - Q{k - A)] , (A7) 

and so for ip. Here is the unit step function. Equation ([6]) reads 

ScTp[(p, v] = ScTp[(p<, V<] + A'S'[0<, v2<, (f>>, V5>]. (A8) 

After integrating out the modes within the shell, the effective action for the surviving modes is 

5'cTP<[0<,V5<] = 5'cTp[0<,</'<] +55'[0<,v9<], (A9) 

where 

e^^^[*<'^<l = j V^^V<P^ exp (^2A^[0<,y.<,0>,y.>]^ p[0>(k, 0), y.>(k, 0)]. (AlO) 
This is essentially Eq. The task here is to compute 5S to order 5s. 
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1. Computing SS to order 6s, perturbative approach and Feynman rules 

We will compute 6S perturbatively. To do this, we will separate in AS", Eq. ( lASp . a free 
CTP action 5*0, which will define the propagators to be used in the Feynman diagrams, 

AS[(f)<, v?<, 0>, V5>] = So[(f)>, (p>] + Si[4><, v9<, 0>, v3>]. (All) 

Momentarily and for sake of brevity, we will write 0> meaning the pair {0>, and (/)>o 

instead of 0>(k, 0) (analogous meaning for 0<). Equation ( lAlOp reads 



e'''^^<^ = j V<p^ expt\^So[<P>] + ^/[0<,0>]J p[0>o]. (A12) 
This is identical to the usual expression for a vacuum-vacuum amplitude. Hence, 

j P0> e^^°[<^>] p[0>o] j - ^ E ^ [<^<] (A13) 



where 

(^?)c [0<] 



I?0> e^^«[^>] p[0>o] 5/[ 



I?0> e^^°[<^>] p[0>o] 



(A14) 



connected 



is the sum of the connected diagrams drawn from 5*7, with the fields 0< acting as external 
currents. The first, constant term in (1A13P can be omitted. Thus, essentially, we can write 
(restoring the full notation for and (f) 

5S[0<, y.<] = {Sj)^ [0<, v.<] + ^ {S])^ [0<, ¥.<] + ... (A15) 

Internal lines in the diagrams will correspond to the propagators deduced from So[(j),(p]. There 
are two propagators 

(V9(k,t)0(k',t'))o = -tGik,t,t')6\k + k'), (A16) 

{^(k, t) y.(k', t'))o = Gik, t, t') + k'). (A17) 

The subscript means that the expectation values are taken with respect to the free action 5*0. 
In general, the function G is zero if t < t'. The third propagator {4>(j)') is always zero. These 
are two general properties of the CTP formalism. 

It is clear that diagrams with more than one loop can be ignored to order Ss. This is because 
for each loop there would be one independent momentum integration over a region of volume 
of order Ss. Thus, the resulting term would be (at least) proportional to {Ss)^, where L is the 
number of loops in the diagram. Hence, it has to be L < 1. In fact, diagrams with L = and 
L = 1 are both at least of order Ss. Owing to the fact that internal lines carry momenta in 
the infinitesimal momentum shell, when computing SS to order Ss, it is enough to consider two 
types of diagrams 0] , shown in Fig. [2] and Fig. [31 
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FIG. 2: Type 1 diagrams have one loop with momentum fc = A. The total external momentum at each vertex 
is zero. 




FIG. 3: Type 2 diagrams are tree diagrams with no external momentum entering at the intermediate vertices. 
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Diagrams of the first type are shown in Fig. [2l They are one-loop diagrams. External lines 
attached to each vertex must have total momentum equal to zero: for each vertex with external 
lines carrying momenta qi, q2, . . . , there will be a Kronecker delta 

^Q;0, (A18) 

where 

Q = qi + q2 + . . . (A19) 

Thus, the same momentum k = AQ, with \Q\ = 1, traverses the loop at every point. The 
modulus of the momentum in the loop is fixed, but not its direction Cl. Hence, instead of an 
volume integration over k, there will be just angular average over Cl and a multiplicative factor 
accounting for the volume of the shell, that is 

n.A'Ss J ^. (A20) 

The fact that external lines at each vertex must carry total momentum equal to zero can be 
understood graphically. Consider for example the one-loop diagram in Fig. 4. 




q2 q4 



FIG. 4. One-loop diagram used to illustrate the fact that external lines should carry zero total momentum. Here 
Q = qi + q2. Both k and k + Q have to be in the momentum shell. 

The analytic expression for this diagram will include an integration of the form 

J cilGi(A:)Gi(|k+Q|)/(qi,q2,q3,q4,k,-k-Q), (A21) 

where Q = qi + q2 is the external momentum entering at the left vertex, and Gi and G2 are 
propagators for modes of the field in the momentum shell bA < k < A. Then, two conditions 
will determine the integration volume, namely 

bA<k<A (A22) 

and 

6A<|k+Q|<A. (A23) 
These conditions are depicted in Fig. [51 For 6s <^ 1, the volume in which the two conditions 
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FIG. 5: In d = 2, the integration region for the diagram in Fig. 4 is the intersection of the two shells. Its volume 
is of order Ss'^ unless Q ~ Ss. 

are satisfied is of order 5s^, except wlien Q is itself of order 6s or superior. In this case the two 
shells coincide with a precision at least of order 6s, and the integration volume is of order 6s. 
In the limit in which 6s 0, the diagram will be null unless Q = 0. Then, if k = fcfi, we can 
rewrite Eq. flA2ip as 

/(qi, -qi, qs, -qs, A^^, -A^^) Ss. (A24) 

Diagrams with more than two vertices can be analyzed in similar terms. 

There are two other rules for the diagrams of the first type. For each internal line connecting 
fields at times t and t', there will be a factor Q{A,t,t') for contractions of the type {(p^p), or 
—iG{A,t,t') in the case of {(p4>). For each internal line there will be two time integrals between 
and T, associated with the times of the two fields of type > joined by the line. 

Diagrams of the second type, shown in Fig. [3], consist of a single chain of propagators, all 
carrying the same momentum k, with |k| = A. For each internal line connecting fields at times 
t and t', there will be a factor Q{A,t,t') or G(A,t,t'). For each intermediate vertex (i.e., not at 
the extremes of the chain) with external lines carrying total momentum equal to q, there will be 
a Kronecker delta 6q-o (the argument is essentially the same that was given above for diagrams 
of the first type). Finally, if the total external momentum entering to the diagram at one of its 
extremes is k, there will be a factor A6{k — A). As before, there will be two time integrals for 
each internal line. 

Note that in both types of diagrams, two internal lines at most are attached to each vertex. 
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APPENDIX B: PROPAGATORS 



The propagators obtained for the free action Sq of Eq.( |T8|) . with initial conditions given by 
Eq. (1251) . are defined as in Eqs. (1A16P and (]A17p . with 



G{k, t, t') = - sin [uj{t - t')] e-"(*-*') e(t - t'), (Bl) 

OJ 



and 

Q{k,t,t' 

2 



a 



with 



and 



[1 + 2/(a /?)] - 7^ U 4 cos[u;(t - t')] - 4 cos[c^(t + t')] + - sm[uj{t + t')] ] 
+ ^ ( \cos[Lj{t - t')] + - sm[u;{t - t')]] e-'^^*"*') Q{t - t') + {t ^ t')\ 
+ Dl [1 + Via 13)] {cos[cu(t - t')] - cos[a;(t + t')]} (B2) 

= A;2 + m^ (B3) 
= ujI~k^, (B4) 

/(x) = (e^' - 1)-^ (B5) 



APPENDIX C: OBTAINING 5S 



Here we will sketch how to compute the actual values of r], (5k, etc. from Eq. ( ]A15I) . We 



will follow the procedure presented in Appendix |Al To illustrate the method, we will give 
some examples of diagram calculations, and show how to extract from 6S the corrections to the 
parameters in the free action and the value of rj. 

According to the definitions ( 1A8I) and ( fTTl) . and due to the quadratic nature of S'o, it results 

AS'[0<, V9<, 0>, V2>] = 5'o[0>, ip>\ + 5'int[0< + 0>, <^< + - 5'int[0<, ^<]- (CI) 

Then, from definition ( lAlip . it is 

5'/[0<, V5<, 0>, V5>] = 5'int[0< + 0>, V5< + V5>] - 5'int[0<, V5<]. (C2) 

We will take as an example the term vn (pip^. Writing = 0< + 0> and = y9< + it yields 

4>l<f2<f3f4: = 4>l>'P2><f3>f4> + 3 (j)i>(f2>f3>f4< + 4>l<f2>f3><f4> + 3 (f)i>lf2>f3<f4:< 

+3 4>l<f2<<-P3>f4.> + 4>l>f2<f3<f4< + 3 (j)l<<-P2<f3<<-P4> + 4>l<f2<f3<f4<- (C3) 
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Here 0i< stands for 0<(ki,ti), and so on. The fact that the function vn is symmetric with 
respect to the variables corresponding to the fields accounts for the factors 3. To write 5*/, the 
last term in Eq. flC3l) . which corresponds to the term S'int[0<, V'<] in Eq. flC2l) . must be taken 
apart. Moreover, since, as it was said before, at most two internal lines can be attached to any 
vertex, we can discard the terms in Eq. ( 1C3I) with more than two fields of the type >. Hence, 
corresponding to (jy^fi^ in 5'int, Sj will display the following terms 

^"d" {\\ j dti j (fk^j 5"* (ki H h k4) 0i>v?2>V53<V54< + 3 (l)i<(f>2<Vi>^i> 

+01>V52<V23<V54< + 3 01<V52<V53<<^4>)f4l({k}; {t}) . (C4) 

Using solely these terms, we can construct two diagrams with one vertex (shown in Fig. 6) and 
which correspond to the term (S*/) in Eq. (1A15I) . and four diagrams with two vertices (shown 
in Fig. 7) corresponding to the term [S]) in Eq. (1A15I) . 

Since all the propagators will be evaluated at A, we will write G(t,t') for G(A,t, t'), and so 
for Q. 





Q 

SI 






-iG 

SI 




0< 


Diagram l.I 






Diagram l.II 





FIG. 6. One-vertex diagrams constructed using the first two terms in Eq. (jC4[) . Continuous (dashed) lines 
represent ip (</>) fields. 



32 



0<1 






0<1 


0<2 




-iG / 




\ 


^ / 








\ 










\ 


\ > 








\ . 


' \ / 










V^^41 




/ ^ \ 






c \ 




1 liQrrT*QTvi / 1 
J-Zldgl dlli Zt.i 




Dine 


!,! dlli Zt.ii 








0<i 


V'<4 


\ ~^ / 


\ 










v 


/ 








^^4l\ 


Ml 










4— <^<5 










-iG \ 




/ fiY"" \ 








V5<2 








V5<6 




Diagram 2. Ill 




Diagram 2. IV 



FIG. 7. Two- vertex diagrams constructed using the last two terms in Eq. (jC4p . Continuous (dashed) Unes 
represent {(f>) fields. 

We use continuous lines for (f fields and dashed for 0. In this way, internal lines corresponding 
to propagators Q are represented by continuous lines, and corresponding to propagators —iG 
by half dashed, half continuous lines. We use 0<j for 0(kj,ti) and so on. Now we analyze each 
diagram in separated subsections. 



1. Diagram l.I 

The first one- vertex diagram in Fig. 6 comes from the term 3(/)<(/9<(^>(^> in Eq. flC4p . after 
contracting the last two fields. There just one way in which this contraction can be made, 
so there is no additional combinatorial factors. According to the rules enunciated at the end of 
Sec. lA 11 there would be a Kronecker delta 5|ki+k2|;05 but the condition |ki + k2| = is already 
granted by the Dirac delta of momentum conservation. Thus, in this case the Kronecker delta 
can be omitted. The contribution of this diagram to 6S reads 

[ dti [ dt2 [ d'^k(f){k,ti)yD{-k,t2) v{k;ti,t2)Ss, (C5) 

Jo Jo JbA 
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where 



v{k;ti,t2 



Jo 



dh / dt4:Vii{\<i,-\<i,tl,-Vt]ti,t2,h,U)Q{h,U) 



We have written v6s and not directly something hke 5f2i'^\ because this is not yet the actual 
contribution of Diagram l.I to 5v2i ■ We have to remove from v those terms associated with 77, 
5m? and 5n. We must separate in v{k] ti, something proportional to 5{ti — ^2) and something 
proportional to d5{ti — 12) /t2- Two projectors are introduced. Given a function of two times 
v{k; , ti, ^2)5 we define 



(C6) 



Vv{k;tiM) = Vv{k) 6{ti-t2) 
and, if v{k; ti, 12) = for ^2 > ^1, 



Qv{k;t,,t2) = Qv{k) 



where 



and 



Vv{k) 



1 
T 



2 ( ^ + 5{t2) - 5(0) ) 5{h - t2) 



dti I dt2 v{k;ti,t2), 
Jo 



Qv{k) 



T 



dti I dt2 v{k;ti,t2) {t2 - h). 
Jo 



(C7) 



(C8) 



(C9) 



(CIO) 



It is easy to verify that = P, = Q, and that QP = PQ = 0. This proves that the 
decomposition 



v{k; ti, ta) = Pv{k; h, ^2) + Qv{k; t^, ^2) + ^v{k; ti,t2 



is unique. Defining 



and 



and using (IClip it yields 



vo = Vv, 



vi = Qv, 



(Cll) 



(C12) 



(C13) 



dti / dt2 / d% (I){k,ti)ip{-k,t2)v{k;ti,t2) 

'0 JbA 



I dti ! d%^(k,ti)ip{-k,ti)voik)+ [ dti [ dt2 [ d'^k (j){k,ti)(p{-k,ti) vi{k) 

Jo JbA Jo Jo JbA 

+ [ dti [ dt2 [ d''k(P(k,ti)y^{-k,t2)Av{k;ti,t2). (C14) 

Jo Jo JbA 
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[Observe that the extra delta factors in (ICSP get rid of the boundary terms, and that the factor 
2 appears because the deltas are evaluated at the extreme of integration.] In this way, we have 
extracted from v{k]ti,t2) two quantities: VQ{k) which acts as a momentum dependent mass 
squared term, and vi{k) which is equivalent to a damping constant. Indeed, comparing with 
(EID, we get 

6^^^-^^ = -vi{k) 5s. (C15) 
From VQ{k) we can extract contributions to and t]. Writing 

„„W=.„(0) + *M5)+i!^ + A„„„), (C16) 



matching terms with Eq. fl27j) . we identify 

Sm^(^-^^ = -2vo{0) 6s (C17) 

and 

After subtracting from v{k; ^1,^2) all the terms belonging to the free action, the net contribution 
from Diagram (l.I) to Sv2i{k;ti,t2) comes from Eqs. (]C11|) and (]C16|) . and is 



6viY\k;tut2) = Av{k;ti,t2) + Avo{k)6{ti-t2). (C19) 

If v{k;ti,t2) were local in time but not constant, that is, if it were v{k;ti)6{ti — ^2), it could 
be interpreted as a time dependent mass. Our definition in this case would give for Sm^^^'^^ the 
time average of f (0; ti) 5s. The definition of P can be generalized for functions g{k; t) depending 
on just one time, writing 

Pg{k;t) = Vg{k), (C20) 

with 

Vg{k) = ^ j\tg{k-t). (C21) 
Something analogous can be said about corrections of the form 

if dti [ dt2 [ d'^k^(k,ti)(f){-k,t2) w{k;ti,t2)Ss. (C22) 

Jo Jo JbA 

We have to subtract from w 6s the part that corresponds to the term 

Si^ (j){k,t)(f){-k,t) (C23) 

in fl27|) . In this case is 

5iy = 2wo{k)Ss (C24) 

where wq = Vw. 
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2. Diagram l.II 

Diagram l.II in Fig. 6 comes from the term 30>V5>V2<V^< in Eq. ( ]C4I) after contracting the 
first two fields. As could be suspected, it is zero. Terms with no fields are not included in the 
definition of the action [Eq. (fTTll and below], and owing to the causal nature of the propagator 
Eq. (]A17ll . they cannot be generated . Remember that we have defined the action in Sec. 
mil in such a way that the there is always at least one field (f) (let us call it cf)*) evaluated at 
a time t* greater or equal than the times of all the Lp fields. When any of the <y9(k, t) fields is 
paired with 0*, there will be a Q{t — t*)^ that means that t must be greater than t* . The two 
conditions are mutually exclusive and the resulting diagram is effectively zero. In the case of 
the vertex v^i there is just one 0, which must be 0* necessarily. 

This is essentially the mechanism which made all the diagrams with no external fields equal 
to zero. 



3. Diagram 2.1 



The first diagram in Fig. 7 contributes to 5f 41 . It has a combinatorial factor equal to 3 (ways 
of choosing which one of the if in the left vertex is x 3 (ways of choosing the in the 
right vertex) x 2 (ways of pairing the fields (v3>, V9>) on the left with the fields (0>, on the 
right) X 2 (combinatorial factor from the expansion of S'j). There is also a factor —i^/2\, which 
comes from the i/2 in Eq. ( ]A15[) and the —i in —iG. The result is 



^^if ({k};{t}) 



18 5s 



'|ki+k2|;0 



f dt[... f dt'^ 
Jo Jo 



a7 



U4i(ki,k2,fi, 



ti, t2, ^3, t'^j 



xG{f^, t[) g{f^, t'^) V4i{n, -n, ks, k4; t[, 4, ts, ^4) 



(C25) 



J 234 



where the subscript ij . . . means symmetrization with respect to the given variables (remember 
that according to our definitions, f„m is symmetrical with respect to the permutations of the m 
fields of type 0, and of the n — m of type ip). 



4. Diagram 2. II 



The second diagram in Fig. 7 contributes to 5f42, has a combinatorial factor equal to 3 x 3 x 2 
and it gives 



S\ki+ka\;o I dt[... I dt'^ I ^V4i(ki,k3,fi, -f];ti,t3,ti,4) 



'y4i(k2, k4, ^2? ^4; '^35 '^4, 



(C26) 



J 12 34 
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5. Diagram 2. Ill 

The third diagram in Fig. 7 is zero due to causahty, which prevents diagrams with no external 
fields. 



6. Diagram 2. IV 

The last diagram in Fig. 7 is an example of a tree diagram. It contributes to Svqi, is 
proportional to 5(|ki + k2 + kal — 1), and has a combinatorial factor equal to 3, which is the 
number of ways of choosing among the 3 ip fields [see Eq. f lC4|) ]. It results 



'54f''H{k};W) = 3 



5(|ki + ks + kgl-l) / dt[ dt'^V4i{ki,k2,k3,-ki-k2-k3]ti,t2,t3,Q 
Jo Jo 



xG{t'^,t[)vii{ki + k2 + k3, ki,k5,kQ]t[,t2,t3,t4) 



(C27) 



23456 



APPENDIX D: RESCALING 



To see in detail the effect of the rescaling, take a generic term Vnm in the interaction part 
(terms in the free action can be analyzed in the same way). After step i) has been performed 
the corresponding term in the action is (constant factors omitted) 

[ dti... [ dtn ! d% ...I d% (5'^(ki + ■ ■ ■ + k„) 

Jo Jo Jbk JbA 

Xv'^^iki, . . . ,kn]ti, . . . , tn) [0(ki, ti) ... 0(km, tra)] [(/^(km+l, Wl) • • • V(kn, ^n)] , (Dl) 

where v'^^ is given in Eq. fl28l) . Then, redefine the fields and change integration variables 
according to Eqs. ( l29|) -(!32l). The momentum integrals are again restricted to < A = 1, but 
the time interval goes up to 6"* T: 



j dti... dtn j d% ...I d% 5\ki + ■ ■ ■ + k„) 

Jo Jo J A J A 



x<m(&ki, . . . ; 6 . . . ) [0(ki, ti) . . .] [<^(k„,+i, t 



m+l ) 



(D2) 



We expand f^^(6ki, . . . , 6k„; 6 "*t„) around kj and ti, and the time integrals as 

functions of their superior extreme around T. To order 6s we get 



dti... dtn d%... d'^kn S'^iki + ■ ■ ■ + k„) 

Jo J A J A 

d d \ 

1 + <j -d{n - 1) + nat - na^ - kj^ + atti — \ 6s 



dti 



^rim(kl5 • • • 1 ^1) 



X [0(ki,ti)...] [^{k 
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-atT 5s y2 dti... dtn I d% ... i d% S'^i^i + ■ ■ ■ + k, 

. ^ Jo Jo J A J A 



x26{tj-T) v'„Jk,,...-t,, 



(ki,ti: 



m+l 5 



(D3) 



Because of the 6{tj — T), the last two hnes are a sum of terms with n — 1 time integrals. One 
field (the one with index j) in each term of the sum is evaluated at time T. Thus, in principle 
boundary terms could appear in the action. But this is not the case: if j < m, the field evaluated 
at tj = T is of type 0, and the corresponding term is zero, because 0(k, T) = by the CTP 
condition ( fT6|) . If j > m then the field evaluated at the boundary if of type ip. But we have 
imposed the condition that in each term of the action there is always a field evaluated at a 
time t* which is equal or greater than the times at which are evaluated all the ip fields. Because 
tj = T, t* must be equal to T. But, again, the CTP condition implies that the term will be zero. 
Thus, we could rewrite Eq. (ID2p keeping T as the superior limit of the time integrals, which 
means that T can be considered as an invariant quantity with respect to the RG transformation. 

In conclusion, rescaling takes the resulting parameters after mode elimination, Eqs. (127|) and 
( l28ll . and returns (to order 5s) 



m? + 5m? + ~ 2a(^ + at — djm^ 5s, 



(D4) 



— — 4)-^' 

u + 5^ + [ -2a^ + at - d - k— ) u 5s, 



7/ -,\ 1 ^ 9 
-na^ + nat-d{n-l) - + atti — 



Rescaling also affects the density matrix. After rescaling it results 
p'[0(k,O),v9(k,O)] =exp 



(D5) 
(D6) 
(D7) 



4 



tanh I "'"^f"*') ^(k,0)^( 



-k,0) 



+ coth(«^l) ^(k.O)^(-k.O) 



exp <- d k 



A 4 



, \a(bk)b^+^"^]\mk)b-'^-^"^] , , , , 
tanh I ^^-^^ i l(/)(k,O)0(-k,O) +. 



■. fD8) 



Then, rescaling takes a and /3, and returns 



d 

a + ( —d — 2a^ — k— ] a 5s, 
f3+ { d + 2a^-k-^j (3 5s. 



(D9) 
(DIO) 
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1. Final remark 



Before leaving this section, we note a special case of rescaling. In principle, one does not 
have to rescale the arguments of (5k, and 5vnmi because they are already of order 5s. Note 
that there are no fc-derivatives nor t-derivatives of these quantities in the differential equations 
of Sec. IIVCI But it may happen (for example in the case of tree diagrams) that 

5t;„„({k};{t})oc5(|K|-A)5s, (Dll) 

where K is a linear combination of the kj. Rescaling kj gives 

5vnm{b{\i}; {t}) oc 5{b\K\ - A) 5s (D12) 



or 



5t;„^(6{k}, {t}) oc (5(|K| - 6-1) 5s. 



(D13) 



The first h ^ can be replaced by 1, but h ^ ^ 1 + 5s inside the delta must be conserved because 
it can define how to take limit values of the functions. We define 5{k — l"*") such that 



r dk5{k-K^ 
Jo 



0, if g < A, 
1 if g > A 



(D14) 



and then rewrite Eq. flD13P as 

5t;„^(6{k},{t})oc5(|K|-A+) 5s. 



(D15) 



APPENDIX E: CLOSED SET OF COUPLINGS TO ORDER 



Here we find the set of couplings for which the RG transformation is closed to order A^. 
Starting from the usual theory, Eq. (!22l) . new terms generated by the RG, and v^i and 
themselves, will not longer be local in time. After one infinitesimal step, there will be terms 
with fields evaluated at two times; for example the one generated by the diagram in Fig. 8. 




FIG. 8. One-loop diagram drawn using two identical vertices local in time, with coupling constant A. The 
resulting term in the effective action will be no longer local in time, but will couple fields at t {4>ip) and t' {ip'ip'). 
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We used and ip for fields evaluated at time t, and 0' and </?' at t' . If we iterate the transformation 
once more, we will get terms with fields evaluated at 3 and 4 times. We may suppose that, 
basically, V41 and W43 remain local and time independent, and that the non local and time 
dependent terms are corrections. 

We assume that ^41 and v^^ are of order A and compute the RG equations to order A^. We 
have to find what kind of couplings can be generated when proceeding up to this order. We 
assume that for s > in the interaction part of the action there will be two terms local in time, 

I dt I d'%i...d'%4 5'^(ki + --- + k4) y4i(s)0(ki,i)^(k2,t)(/?(k3,i)(^(k4,i) (El) 

and 

I dt I d% . . . d% (ki + • • • + k4) ^43(5) 0(ki, t)0(k2, t)0(k3, t)<^(k4, t), (E2) 
Jo JbA 

where V41 and V43 are order A and do not depend on {k}, and satisfy 

y„(0) = y,3(0) = -j|^. (E3) 

Starting from these vertices we perform successive infinitesimal RG transformations until no new 
terms are generated. The finite RG transformation will be closed with respect to the resulting 
set of couplings, a desired property which was stressed before. We do not need to keep trace of 
the precise value of each coupling from step to step. At this point we just want to enumerate 
the couphngs that have to be considered to be consistent to order A^, and find their general 
functional form. Rescaling have to be taken into account only in the deltas appearing in tree 
diagrams. 

Then, for the first step, wc start with V41 and V43 and construct all the possible diagrams 
of order A and A^, shown in Fig. 9. They all will give new terms. Below each diagram we 
have written schematically the corresponding term, indicating the name of the coupling and the 
variables on which depends. We have defined 

<5i2... = |ki + k2 + ...|. (E4) 

The couplings have, by definition, the same symmetry that the set of fields that they couple. 
Most diagrams give a result that has already the required symmetry. 

Next, for the second step, we combine the original terms with the new ones and construct 
all possible diagrams of order A and A^, shown in Fig. 10. Only two new terms are generated 
at this step, with couplings W21 and W22, all of order A^. They cannot be used to produce new 
terms, because combined with any other vertex the resulting diagram would be at least of order 
A'^. Hence, a third step fails to create new terms. In Fig. 10, we have drawn the diagrams for 
the two field couplings at the end, because, in some sense, they are at the deepest level and 
depend on all the terms which precede them. Note that diagram (8) in Fig. 10 will not depend 
on k, owing to the fact that W41 is evaluated at Q = 0, and hence the dependence on k cancels 
out. 
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FIG. 9. Set of diagrams obtained after performing one infinitesimal step of RG transformation starting from 
the initial condition given by Eq. (|22p . Below each diagram we indicate the coupling generated, with its 
corresponding dependence on the momentum and time variables, with Qij,,, = |ki + k2 + . . . |. 
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FIG. 10. Set of diagrams obtained after performing two infinitesimal steps of RG transformation starting from 
the initial condition given by Eq. (j22p . A third step does not create new couplings. 
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Some conclusions can be drawn. If m^, n and u are initially zero, will be of order A through 
V41 in diagram 1 of Fig. 9, and k and u will be at least of order A^, through the one vertex 
diagrams of Fig. 10. Note, as well, that rj will be of order A^ owing to the fact that 14^21 is order 
A^ and is the only term of the type (f)(p which depends on k. This means that, at order A^, rj can 
be omitted and at can be set equal to 1 in the left-hand side members of Eq. (!38l) . and Eqs. 
(l35l) - (!371) . We will assume that this is the case. 

The vertices which depend on two times, have half of the fields evaluated a t and half at 
t'. Diagrams which depend on two times with an uneven number of fields at t and t' can be 
drawn using a fgj vertex, but they are zero. This fact can be understood noting that, after 
iterations of the RG transformation, VQi{k; t, t') will be given by a sum of terms, with the nth 
term proportional to 5{k — (remember that A = 1). So, if k is not in the interval &~^], 
VQi{k]t,t') is zero, or in other words (and it will shown below) 

vei{k; t, t') oc Q{k - l+)e(en+ - k), (E5) 

where we use l"*" meaning that Q{k — 1+) is strictly zero if A; = 1. Then, consider for example 
the case of 

V6l(|ki + ks + k3|;t,t') Lf' Lf' if' . (E6) 

If we connect the two fields to generate a term of the form 

y^W, (E7) 

we must set ^^2 = — ka = 1 and integrate over Vt. But because of (lESp . the unpaired field 
should have momentum 



|ki| > 1, (E8) 

which is impossible, because ki < 1. The same occurs if we connect two of the if'. The 
conservation delta, (5'^(ki + k2 . . . ), which is always present, allows us to write the first condition 
on Eq. (IE5]) as 

e(|k4 + k5 + k6|-l + ). (E9) 

Setting, for example, ks = —kg = 1, it should be |k4| > 1, and the conclusion follows. 

There are other possible couplings depending on t and t', with an even number of fields at 
each time, that due to causality are not generated to order A^. For example, from 

t;g)000y.VV (ElO) 
joining one of the (f) with one of the y?' we cannot generate a new term of the form 



00 v'V, 



(Ell) 
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because, on one hand, necessarily is f <t (see Sec. IIIip . and on the other hand, the propagator 
G{t',t) connecting cp' with entails the opposite condition t' > t. Also, vq^ (fxjxf) 0VV' cannot 
create a term of the form (fxp (f)'(j)', because by construction (but not necessarily, as in the previous 
example) is t' <t. 

Other typical cancellation occurs in the diagram of Fig. 11 

'^^ oce(t-t') 

V^43 4 fV21 =0 

oce(t'-t) 

FIG. 11. Example of a diagram which is null due to a closed chain of steps functions of the form: Q{ti —^2)0(^2 — 
ts) . . . Q{tn-i — tn)Q{tn — ti). This is typical of diagrams drawn using vertices which are local in time. 

There are other cases that can be analyzed in similar terms. 



1. Final remark 

We have seen in Sec. Othe method used to extract from diagrams with two external lines 
the terms which by definition belong to the free action, Eq. (fTSl) . Something analogous must 
be done with diagrams which give terms with 4 external fields of the form 

(E12) 

and 

(j)ip(f)'(f)'. (E13) 
For example, the contributions of diagrams (2) and (5) of Fig. 10 can be written as 

[ dt [ dt'l d%...d%5^{ki + ...) </.(ki,t)(^(k2,t) ip{ks,t')ip{ki,t') v{Qu;t,t')Ss, (E14) 

Jo Jo JbA 

f dt f dt' f d%...d%S'^{ki + ...) 0(ki,t)(/?(k2,t) 0(k3,t')0(k4,O w{Qi2;t,t')Ss, (E15) 

Jo Jo JbA 

where 

v{k- 1, t') = 18 6k-o Vl Gg{t, t'), (E16) 

and 

w{k;t,t') = 18<5fc;o V^iV^s Gg{t,t'). (E17) 
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In principle, v 5s and w Ss are not directly identified with SW^i and but can contain 

contributions to 6V41 and 6V43. Using the definition (1C9I) . we write 

v{k; t, t') 5s = 6{t - t') 5V41 + 5W4i{k- 1, t'), (E18) 

w{k; t, t') 6s = 6{t - t') + dW^sik; t, t') (E19) 

with 

SVii = Vv{0) 6s, (E20) 

and so for V43. In this way, the corrections to V41 and V43 are isolated, while 6W41 and 5^43 
give the net contributions to the W41 and W43. 

APPENDIX F: DETAILED EXPRESSIONS FOR k AND u 

Here, the formulas for k and u, with = are given. The expressions are for k = 0, but 
are asymptotically valid for < k when 1 z. They are composed by two parts, corresponding 
to 2; < 2 and z > 2, which join smoothly for z ^ 2. 

k{z, T) = K<(2, T) 6(2 -z) + K^{z, T) Q{z - 2), (Fl) 
u{z,T) = u^{z,T) e{2 - z) + u>{z,T) 9(^-2). (F2) 



1. Expressions for k 



We write 
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«^< = K)'5Z«^<i, (F3) 

i=l 

and an analogous equation for where 

= Aj _63z log3- -(3 -4^ + 2^) + 6(-5 + 3;z)cosr- 14(-1 + z)cos(3r)|, (F4) 

18 r 2 

K<2 = — I - COS [r(-2 + ^)] + - cos(T2) - 7 cos [T(2 + ^)] 

+ cos[T(l -2^)] + 7cos[r(l + 22)] |, (F5) 
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9 r 

^c^^= )z log(-2 + ^) - 22;(5 + 4 cos T) log z-7z log(2 + z) 

-4zlog(-l + 2z) + 28zlog(l + 2z)}, (F6) 

«;<4 = -18 [z sin T - sm(rz)] , (F7) 

«;<5 = 3| _ 3 (1 + 3^) Si(r) - 21(-1 + Si(3r) - 3(-2 + Si[r(-2 + z)] 
+12 Si(Tz) - 21(2 + z) Si[T(2 + z)] + 3(-l + 2z) Si[T(-l + 2z)] 

+21(l + 2z)Si[T(l + 2^)]}, (F8) 



«;<6 = 18^T [Ci(T) - Ci{Tz)] + ^{ - 3Ci(T) - Ci[T(-2 + z)] + 7Ci[T(2 + z)] 

+21Ci(3r)+4Ci[r(-l + 2z)] -28Ci[r(l + 2^)]}, (F9) 



= —j - 12[Ci(2r) - Ci(2rz)] COST + 3 [Ci(r) - Ci{Tz)] cos(2r) 

16 [Si(2T) - Si(2T^)] sin T - 4 [Si(T) - Si(T^)] sin(2T) | , (FIO) 

201 ^63 2,z ( „ 
'^---8T-^^+8T;+32Tr 
-24[7£ + 211og3 + log64 + logT- (2 + log4)T2]}, (Fll) 

= TTT^^f-z^ [l44cosr-24(l + log4) cos(2T) - 112cos(3T) - 105cos(4T) 
32r2; y L 

-28cos(2rz) + ^| - 240cosr + 120cos(2r) + 112cos(3r) - 84cos(4T) 

-24 cos[2T(-l + z)] + 48 cos[T(l - 2z)] + 336 cos[T(l + 2^)] }^ , (F12) 

/^>3 = -^|^[8sinr - 4sin(2T) - 7sin(4T)] + 14sin(2T^)|, (F13) 
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9z 
"4T 



9z 

«;>4 = - — [2(9 + T^) + 8 COST - cos(2T)] \ogz 
- log(-l + z)+A log(-l + 2z) - 28 log(l + 2z)] , (F14) 



9z r 

«;>5 = 4j^{(-3 + 4T')Ci(T) + Ci(2T) + 21Ci(3T) + 7Ci(4r) - Ci[2T(-l + z)] 
-2T^ [2Ci(2T) + 7Ci(4T) - 701(27^)] + 4Ci[T(-l + 2z)] - 28Ci[T(l + 2z)]}, (F15) 

= ^1 _ (1 + 3^)Si(T) + (3 + z)Si{2T) - 7(-l + z)Si{3T) - 7(1 + z)Si(4T) 
+(1 - ^)Si[2r(-l + z)] - 7Si{2Tz) + (-1 + 2^)Si[r(-l + 2z)] 

+7(l + 2^)Si[T(l + 2z)]}, (F16) 

9z r r 1 
ft;>7 = — <^ 12Ci(T) - 12Ci(2r) - 2Ci(4T) + 2Ci{2Tz) cos(2T) 
ST L L J 

+64 [Si(2r) - Si{2Tz)] sinT - 2 [8Si(r) - 8Si(2r) + Si(4r) - Si{2Tz)] sm{2T) 

+48[-Ci(2r) + Ci(2r^)]cosr|. (F17) 

The functions Si and Ci are the sin and cos integral functions, and je is the Euler's constant. 

2. Expressions for u 

We write 

7 



i=l 

and an analogous equation for i/>, where 



iy<i = ^ [-12 - log 9 - 32z + (44 - log 2187) z^]), (F19) 



i/<2 = 27 z [-Z sin T + sm{Tz)] , (F20) 

i/<3 = ^{ - 3(1 + 7z^) COST - (-1 + z^) cos(3r) + (-2 + z) cos[T(-2 + z)] 
+lScos{Tz) -{2 + z) cos[r(2 + z)] + 7(-l + 2z) cos[r(l - 2z)] 

+(l + 2^)cos[T(l + 2z)]}, (F21) 
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i/<4 = -{ - 3(1 + 9^2) Si(T) - 3(-l + z^) Si(3T) + (-2 + 2)' Si[T(-2 + z)] 
+24z Si{Tz) - (2 + zf Si[T(2 + - 7(1 - 2zY Si[T(l - 2z)] 

+(l + 2^)2Si[r(l + 2^)]}, (F22) 

i/<5 = [-2 + (-19 + 4T^)z^] Ci{T) + (2 + 7z^) Ci(3r) 

+(-4 + z"^) Ci[r(-2 + z)] + 12(2 - T^z^) Ci(r^) + (-4 + z^) Ci[r(2 + z)] 

+14(-1 + 4^2) Ci[r(-1 + 2^)] + 2(1 - 4^2) Ci[r(l + 2z)] }, (F23) 

i^<6 = ^{(4 - log(-4 + z') - 2(12 - 25^2) log^ 

-2(-l + 4z^) [71og(-l + 2z) - log(l + 2^)] }, (F24) 

9^2 r 

= 15777124 [Ci(2T) - Ci(2T^)] COST + [Ci(T) - Ci(T^)] cos(2r) 

+32 [Si(2r) - Si(2Tz)] sin t} , (F25) 

i^>i = { - 4 (7 + 47E + 2 log 1536 + 4 log T) - 56z 

+ [135 + 47i5 + 28 log(4/3) + 4 log T] 1 ^ (p26) 

i/>2 = ^1 - 22sin(2T^) - [12 sin T - 6sin(2T) - sin(4T)] z'^^ (F27) 

1/^3 = -^1 - 24(1 + 7^2) COST + 4(14 - Az + \\z^) cos(2T) 
16T I 

-8(-l + z^) cos(3T) - (8 + 3^2) cos(4T) - 56(-l + z) cos[2T(-l + z)] 
-12 cos(2r^) + 56(-l + 2z) cos[r(l - 2z)\ + 8(1 + 2z) cos[r(l + 2^)]|, (F28) 

= ^{3 [-2 - (19 - 4T2)z2] Ci(T) + 2 [5 - (-7 + eT^)^^] Ci(2T) 
+ (2 + 7^2)01(37) - (2 - ^2 + 4r^z^)Ci(4r) + 14(1 - z^)Ci[2r(-l + z)\ 
-(2 - 4r^z^)Ci(2Tz) - 14(1 - 4^^)Ci[r(-l + 2z)\ + (2 - 8z^)Ci[r(l + 2z)]}, (F29) 
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i/>5 = I - 3(1 + 9z'^)Si{T) + 2(7 - 2z + 7z'^)Si{2T) + 3(1 - z^)Si(3T) 
+ (-4 - ^2)Si(4T) - 14(1 -2z + z^)Si[2T{-l + z)] - 4zSi{2Tz) 
-7(1 -Az + 4z2)Si[r(l - 2z)] + {1 + Az + Az'^)Si[T{l + 2^)] |, (F30) 

i^>6 = + z') log(-l + + (1 + 17z^) log^ 

-(-1 + 4z^) [71og(-l + 2z) - log(l + 2^)] }, (F31) 

i^>7 = ^{24 [Ci(2T) - Ci(2T^)] cosT 
+ [Ci(T) - Ci(2T) - 7Ci(4T) + 7Ci(2T^)] cos(2T) 
+32[Si(2r) -Si(2T^)]smr + 7[-Si(4r) + Si(2T^)]sin(2r)|. (F32) 
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